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Abstract. In the framework of the model that uses black holes endowed with electromagnetic structure (EMBH) as 
the energy source, we study how an elementary spike appears to the detectors. We consider the simplest possible 
case of a pulse produced by a pure e"^e~ pair-electro-magnetic plasma, the PEM pulse, in the absence of any 
baryonic matter. The resulting time profiles show a Fast-Rise-Exponential-Decay shape, followed by a power-law 
tail. This is obtained without any special fitting procedure, but only by fixing the energetics of the process taking 
place in a given EMBH of selected mass, varying in the range from 10 to 10'^ Mq and considering the relativistic 
effects to be expected in an electron-positron plasma gradually reaching transparency. Special attention is given 
to the contributions from all regimes with Lorentz 7 factor varying from 7 = 1 to 7 = 10** in a few hundreds 
of the PEM pulse travel time. Although the main goal of this paper is to obtain the elementary spike intensity 
as a function of the arrival time, and its observed duration, some qualitative considerations are also presented 
regarding the expected spectrum and on its departure from the thermal one. The results of this paper will be 
comparable, when data will become available, with a subfamily of particularly short GRBs not followed by any 
afterglow. They can also be propedeutical to the study of longer bursts in presence of baryonic matter currently 
observed in GRBs. 
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1. Introduction 

It is by now clear that the Gamma-Ray Bursts phenom- 
ena, here generally called GRBs, far from representing a 
single physical event, are quite complex and are compos- 
ite of very different phases, corresponding to distinctively 
different physical processes. Any data analysis procedure 
should be done taking into due account these different 
epochs and avoiding time averaging procedures on differ- 
ent epochs. 

In recent years much attention has been devoted to 
analyzing GRB afterglows, originating from the interac- 
tion of a relativistic expanding fireball with the surround- 
ing baryonic material, leading to important correspon- 
dences between observational data and phcnomenologi- 
cal models 



ing the structure of the burst itself. In the existing litera- 
ture, it has been alternatively assumed that GRBs origi- 
nate: 

a) by the smooth transition to transparency of an opti- 
cally thick electron-positron plasma, whose origin is not 
discussed (see e.g. Goodman, 198(]| ); 

b) from external shock processes in the collision of an ex- 
panding fireball, whose origin is not discussed, with bary- 
onic matter at rest (see e.g. Meszaros & Rees, 1992); 

c) from the interaction with each other of different parts 
of the fireball due to different expansion velocities, by in- 



ternal shock processes (see e.g. Paczyhski & Xu, 1994); 
d) from a mixed "internal-external" scenario, in which 



both effects are present (see e.g. Piran, 199E) 



(see e.g. 



Vietri, 1998i |Sari fc Esin, 2000| and 
Djorgovski, Frail, et al., 2000|) . 

Less progress has been made in understanding the ba- 
sic mechanism which produces the fireball and in explain- 
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e) from a Compton-drag process, caused by the inter- 
action of a relativistic fireball, whose origin is not dis- 
cussed, with a very dense soft photon bath (see e.g. 
phiseUini, et al., 200q ). 

Important as they are, these semi-empirical ap- 
proaches miss the completeness of a comprehensive quan- 
titative model. Some of the above processes could indeed 
be present in the different epochs of the GRBs. In order to 
evaluate their relative magnitude and give a quantitative 
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estimate of the associated rclativistic effects is necessary 
to have a detailed model, with the corresponding equa- 
tions of motion and the time evolution of the system. 

For this reason we consider the model which assumes 
the energy source of GRBs to be the process of vacuum 
polarization around a black hole endowed with electro- 



magnetic structure (EMBH) (see Damour & Ruffini, 1975 
and Ruflini, 1998). For such a system the fundamental en- 



ergetic aspects have been clarified by the introduction of 
the Dyadosphere of an EMBH (Preparata, et al., 1998b) 
as well as the equations of motion for the system have 
been integrated and are summarized in the next section. 

We proceed in four separate steps: 

In this paper, we make a first step. We consider the 
simplest possible case: the one corresponding to a pulse 
produced by a pure e^e~ pair-electro-magnetic plasma, 
the PEM pulse, in the absence of any baryonic matter. 

The equations of motion of this system have already 
been integrated ( Ruffini, et al., 199£ ). We here examine 
how the PEM pulse gradually emits radiation during its 
evolution and the GRB simply occurs during the smooth 
approach to the transparency condition, when the mean 
free path of the photons in the plasma, in the comoving 
frame and A in the laboratory frame, given by 



(1) 



approaches the thickness of the PEM pulse itself. Here 
(j.y_e is the photon and electron Compton scattering cross- 
section, whose Thomson limit is ^^^^ ~ 6.66 • 10~^^ err? 
and rigi is the proper number density of e^e~ pairs. Thus 
the GRB is essentially determined by the electron-positron 
pair annihilation and the expansion and cooling of the 
PEM pulse. 

It has been clear since the classic work of Rees on ex- 
panding radio sources that rclativistic effects are central 
to the u nderstandin g of the astrophysics of extragalactic 
sources ( Rees, 196^ ). The typical Lorentz 7 factor consid- 
ered in that study was of the order of 7 ~ 5 and essentially 
constant over the time of observation of the astrophysical 
source. In the case of GRBs such rclativistic effects are 
also important, but there are three major differences: 

1. The rclativistic effects are much more extreme. It was 



pointed out (see Ruffini, et al., 1999) that, in the case 
of the absence of baryonic matter, a Lorentz 7 fac- 
tor of up to 10^ can be reached. This result has also 
been extended to the case of the expansion of the 
PEM pulse in the presence of baryonic matter (see 
Ruflini, et al., 2000D where even larger values of the 
Lorentz 7 factor can be reached. 

The transition from 7 = 1 to 7 ~ 10^ occurs in a few 
hundreds of the characteristic travel times in the PEM 
pulse (i.e. a few seconds). It is therefore impossible, as 
shown in Fig. ||, in Fig. || and in sections @, ||, |[ @ of 
the present paper, even as a rough approximation, to 
consider the Lorentz 7 factor to be constant during the 
emission process of the GRB. 



3. The transparency condition given by Eq.(|i|) is reached 
gradually by the effect of the cooling and the expansion 
of the plasma and the corresponding decrease in the 
number of the electron-positron pairs. Consequently, 
in principle, all the stages from 7 = 1 to 7 ~ 10* 
contribute to the final burst structure. 

The main point of the present paper is to clarify the in- 
terplay of rclativistic effects at work in this extreme case, 
in order to obtain the elementary spike intensity as a func- 
tion of the arrival time, and its observed duration. In par- 
ticular we here take into account, still for simplicity in the 
case of a spherical geometry, the effects due to: 

— the varying thickness of the emitting region, 

— the energy flux, essentially modulated by a time vary- 
ing "screening factor" (see appendix A and B), 

— the time variation of the 7 factor and consequences for 
the observed arrival time. 

We examine explicitly, for simplicity, the case of Reissner- 
Nordstrom EMBHs with ^ = 0.1 and M = IOMq, 
M = 10^ A/© and M — lO'^M©. Some considerations on 
the expected spectrum and its departure from the thermal 
one are also presented in section (H) and will be further 
examined in forthcoming publications. 

All the results of the treatment presented in this paper 
can be observationally relevant for a very special class of 
short GRBs without any afterglow. They can also be of 
qualitative interest for the first spike and early features of 
a more complex long burst. It is important to stress that 
in the case of the elementary spike here considered the 
entire energy of the Dyadosphere is emitted in the burst. 
This is not the case for the long bursts, the large major- 
ity of the currently observed GRBs, where much of such 
energy is transferred to the kinetic energy of the baryonic 
component. 

In a forthcoming paper (Bianco, Ruffini, Xue, in prepa- 
ration) , we consider the emission from a PEM pulse inter- 
acting with baryonic matter, before the condition of trans- 
parency be reached, and the corresponding consequences 
for the observable effects on the intensity, spectrum and 
time structure of the burst. The equations of motion for 
such a pair-electromagnetic-baryonic pulse (PEMB Pulse) 
have been integrated and the relative intensity to be ex- 
pected for the GRB versus the kinetic energy left in the ac- 
celerated baryonic material (Ruffini, et al., 2000) has also 
been given. 

In a final paper (Bianco, Chardonnet, Fraschetti, 
Ruffini, Xue, in preparation) we analyze the interaction 
of the accelerated baryonic material (ABM Pulse) with 
the interstellar medium. As the transparency condition is 
reached all the baryonic matter is left having acquired an 
enormous Lorentz 7 factor reached at the time of decou- 
pling, typically in the range 10 ^ lO''. The interaction 
of these very high energy baryons with the surrounding 
interstellar medium gives rise to the afterglow. 

We finally consider (Chardonnet & Ruffini, in prepa- 
ration), within the above model, the production of very 
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high energy cosmic rays by the electrostatic acceleration 
process of the remnant EMBH. 



2. Main assumptions of the EMBH model 

The most general black hole expected from gravita- 
tional collapse is an EMBH characterized by a Kerr- 
Newmann geometry endowed with axial symmetry and 
electric and magnetic fields (RufRni & Wheeler, 1971). 
That indeed an EMBH with mass smaller than 7.2 • 
10^ Mq, can give rise to vacuum polarization and to 
the creation of e'^e~ pair-electromagnetic plasma outside 



the horizon via the Heisenber g-Euler-Schwinger proces s 
was clearly demonstrated (see Damour fc RufEni, 1975|) . 
It was there shown how these process can ap- 
proach reversibility, in the sense of Christodoulou and 



Ruffini ( Christodoulou fc Ruffini, 1971 ) and that this phe 



nomenon would lead to a most natural model for GRB. 

The discovery of the afterglow of GRB by the Beppo- 
Sax satellite and the consequent clear determination of the 
distance and energetics of the GRB's has motivated us to 
return to this field. In order to give an estimate of the 
efficiency of the quantum electrodynamical process in the 
gravitational field of an EMBH we have considered the ide- 
alized case of a spherically symmetric Reissner- Nordstrom 
geometry, neglecting all the effects associated with rota- 
tion, to be later examined in a much more complex treat- 
ment. We have then computed the physical parameters, 
the spatial extension, the total energy, and the spectrum 
of the relativistic plasma of electron-positron pairs created 
by the vacuum polarization process and introduced the 
concept of the Dyadosphere (Preparata, et al., 1998a, t). 

The computation of the relativistic hydrodynamics of 
such an e"*" and pairs and electromagnetic radiation 
plasma in the field of an EMBH, still in the simplified 
case of spherical symmetry, have been carried out both 
with simple semi- analytical and numerical approach car- 
ried out in Rome and their validation by the full hy- 
drodynamical numerical computation at Livermore (see 
Ruffini, et al., 1999 ). The main result shows the forma- 
tion of a sharp slab of and pairs and electromag- 
netic radiation: the PEM pulse. Such a PEM pulse keeps a 
constant width in the laboratory frame and in a few hun- 
dreds of characteristic crossing times of the Dyadosphere 
reaches extreme relativistic conditions with characteristic 
Lorentz 7 factors of 10^ - 10"*. The pair density at the 
beginning is so high that the plasma is optically thick, 
and only very few photons can escape from the plasma. 
With the expansion, cooling and annihilation of electron- 
positron pairs, the condition of transparency is gradually 
reached (see Eq.(|^)). The details of the temporal develop- 
ment of the burst emitted as the transparency condition 
is reached, as well some indications on the expected spec- 
tra, are the subject of the present paper. For simplicity, 
we have neglected the feedback of the PEM pulse on the 
adiabaticity condition prior to reaching the final moment 
of transparency. 



The case in which some baryonic matter is engulfed 
by the PEM pulse, prior to reaching the condition of 
transparency, has been the subject of a successive work 
(Ruffini, ct al., 2000). Our simplified model, as well as the 
validation by the full numerical codes at Livermore, have 
shown how the slab approximation is still valid in this 
more general case for a large range of the masses of the 
baryonic matter (a pairs-electromagnetic-baryons pulse, 
PEMB pulse). Most important, the addition of baryonic 
matter leads, through the electronic component, to an 
increase in the opacity of the PEMB pulse. The trans- 
parency condition is now reached at later times, lead- 
ing to an ever increasing transfer of the total energy of 
the PEMB pulse to the kinetic energy of the baryonic 
component. The larger the amount of baryonic matter, 
the smaller is the energy of the PEMB pulse released in 
the GRB, for fixed values of the mass and charge of the 
EMBH. 

All these treatments refer to the idealized case of an 
already formed EMBH, it is clear, however, that in real- 
ity the Dyadosphere will be formed during the process of 
gravitational collapse itself, prior to the formation of the 
EMBH, and such a process may have some distinct de- 
tailed observational signature. The general energetic fea- 
tures of GRBs here presented have been obtained from 
the idealized model which con siders as a starting poin t 
an already formed black hole ( Preparata, et al., 1998b ). 
In order to obtain the fine details observed in the time 
structure of GRBs, as well as their detailed spectral evo- 
lution, there is no alternative but to study the gradual for- 
mation of the Dyadosphere, as the horizon of the EMBH 
is approached and formed. In preparation of this more 
complicated analysis and the one corresponding to axially 
symmetric configurations, we have obtained some prelim- 
inary results concerning the relevant physics of a charged 
collapsing shell in general relativity (Klippert and Ruffini, 
in preparation). The general energetic aspects here con- 
sidered will not be modified but the time constant will 
be longer due to general relativistic effects (Cherubini, 
Jantzen, Ruffini, in preparation). 



3. The radiation flux from the PEM pulse 

The frequency lo and wave-vector k of photons emitted 
from the PEM pulse (see Fig. Q) expressed in the labora- 
tory frame are: 



k = — (— sin -|- cos ^9v) , |k| 
c 



UJ 
C 



(2) 



where d is the angle (in the laboratory frame) between 
the radial expansion velocity and the direction from the 
origin of the PEM pulse to the observer, v is a unit vector 
along the radial expansion velocity of the PEM pulse, and 
u is a unit vector orthogonal to v oriented toward rising 
-d. We are assuming here that k and Rt are parallel, also 
for photons emitted with i? 7^ 0. This is clearly a good 
approximation, because the distance Rt corresponds to a 
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redshift z ^ I, while the radius of the emitting region is 
of the order of magnitude of a few hght seconds. 

Then the Lorentz boost along v to the comoving frame 
of the PEM pulse yields the corresponding comoving quan- 
tities: 



Wo = 70; ( 1 



■ cos-d 



ko = — |k| sin -du + 7 |k| (^cos - 



(3) 



(4) 



Integrating Eq. (|^) over all photon momenta, we obtain 
the total radiation flux in the local laboratory frame of the 
PEM pulse. 

The number of photons (fc —^ k + dk) radiating per unit 
time from a small surface element d^S of the PEM pulse 
to the detector is given by (per unit area of detector): 



'¥k 



{k-^)cd^Y.- 



, dD,]^ 
~A 



d^k = 



A 



(10) 



In the comoving frame photons radiating out of the 
PEM pulse must have (see Eq. 



cos-d > -, 
c 



(5) 



where fc — -r?, is directed towards the observer, A is the 

|k| ' 

area of detector, Rt the distance to the detector from 
the origin of the PEM pulse and dilk is the solid angle 
subtended by the detector. From Eqs.(^,p|), we have 



because the component of the photon momentum in the 
comoving frame along the radial expansion velocity direc- 
tion must be positive in order to escape. 

The large amount of high-energy photon emission is 
mainly due to electron-positron annihilations. Thus, as a 
preliminary consideration and approximation, we assume 
that photons are in equilibrium at the same temperature 
T with electron-positron pairs before and at decoupling. In 
the comoving frame of the photon electron-positron pair 
plasma fluid, the black-body spectrum of photons that are 
in thermal with e"'"e~-pairs is given by 



k -w — cos -d £ 



1 



(11) 



By using Eq.(||), we integrate Eq.([lO|) over all photon 
energies and obtain the observed infinitesimal energy flux 
in the case of an optically thin PEM pulse (in ■^^) ■ 



dF 



T(t) 



'(t) (1 



vit) 



cos-d] 



■c coszJd E. 



(12) 



duy 

IPko 



1 



1 



exp(^)-l' 



(6) 



Integrating only rf0 £ [0, 27r] in the surface-element 
d^E of the PEM pulse, and adding the screening factor 
S (t), we finally obtain the infinitesimal photon energy flux 
directed towards the detector and emitted at time t with 
an angle i? for an optically thick PEM pulse: 



dF {t, I?) = 27ra 



where n-y is the number-density of photons and T is 
the temperature in the comoving frame. The non-thermal 
spectrum due to multiple inverse Compton scattering may 
slightly modify these assumptions. These consideration 
relevant to the details of the spectrum (see section ||) do 
not modify the temporal profile of the total detected radi- 
ation fiux, which largely depends on the screening factor 

S{t) rather than on the radiative spectrum. .^2 ^^-j ccos'dS {t) dcosS (13) 

Note that due to the Liouville theorem on the invari- 
ance of the distribution function (see e.g. Ehlers, 1971), where we omit the factor R^^. dF will be then measured 
we can write: in terms of ^ , and we will still have to multiply the result 

by R^^ to obtain the observed intensities. Note that the 
radiation flux is axially symmetric with respect to the axis 
directed towards the detector. 



Tjt) 

j(t) (1- liMcost?) 
•r^ (t) c cos -d S {t)d COS'S, 



1 



1 



exp 



■. kT 



- 1 



exp 



- 1 



(7) 



where coo and T are the photon energy and temperature 
in the comoving frame, while lu and Tiab are the same 
quantities in the local laboratory frame. We can now use 
Eq.(||) for ujo obtaining: 



Tlah — 



T 



7(1- f cosi?)' 



(8) 



Therefore we can use the following formula for the pho- 
ton spectrum: 



driry 
d^ 



1 



77" / ^0^7(1-^ COS)?) 

exp I — St 



(9) 



4. The arrival time versus emission time 

Due to the high value of the Lorentz 7 factor (lO'^ ^ 10'*) 
for the bulk motion of the expanding PEM pulse, the 
spherical waves emitted from its external surface ap- 
pear extremely distorted to a distant observer (see e.g. 
piees, 1966 ). The surface emitting the photons detected 
at an arrival time difference ta, measured from the arrival 
of the flrst photon, is not trivially the spherical surface of 
the fireball at a fixed time difference t, measured from the 
initial emission, but photons emitted at different t and at 
different angles reach the detector at the same time ta- 
The relation between emission time t and arrival time tp. 
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in the case of a constant 7 for the expanding fireball, has 
been found by Rees (for the definition of ■& see Fig. |l|) (see 
Rees, 1966|): 



tn — t [1 COS 7? 

C 



(14) 



The external radius r (t) of the fireball, due to the con- 
stancy of 7, is obviously given by: 



r {t) = vt 



(15) 



If, in Eq.(n4), we fix the value of the arrival time — 
t*, using Eq.(p^ we can find the equation describing the 
"surface" emitting the photons detected at arrival time t*: 



vtt 



1 - i^cosiJ' 



(16) 



(see 



which describes an ellipsoid of eccentricity 
Rees, 1966| ). 

This is the approximation that has been widely 
used in the gamma ray burst literature: see 



e.g. Sari & Piran, 1997, 


Sumner & Fenimore, 1997 


Ramirez-Ruiz & Fenimore, 1999, Fenimore, et al., 1996 


Granot, et al., 1999 








Fenimore, 1999 


Panaitescu & Meszaros, 199? , 




Piran, 1999 and rcfcr- 



ences therein. 

However, in our case the velocity v is not constant: in a 
few hundreds of the PEM pulse crossing timeQ the 7 factor 
of the bulk motion of the expanding PEM pulse goes from 
to ~ 10'', and so we have to find the generalization of 
Eq.(|l^ for nonconstant velocity. This can be done using 
the geometry of Figure ^ 

We set t = when the fireball starts to expand and 
the first photons are emitted, so that r (0) = rds- Let a 
photon be emitted at time t from the point P. Its distance 
from the observer is L. The time it takes to arrive at the 
detector is, of course, So, its arrival time, measured 
from the arrival of the first photon a time after its 
emission at i = 0, is: 



c c 



(17) 



where we have defined ta = when the first photon emit- 
ted at t = and d = reaches the observer. L is clearly 
given by: 



L = 



i?2 + r {ty - 2 i?T r (t) cos 1} 



(18) 



where at any given value of emission time t, cos")? can as- 
sume any value between and 1 as noted above, where 
V (t) is the expansion speed of the fireball at time t (see 
Eq.(l)). Now, r{t) is of the order of magnitude of some 
light-seconds, and Rt corresponds to a redshift z ~ 1. 



The crossing time of the PEM pulse is given by 



and this is ^ 10 

and Q = O.lQma: 



seconds, for a black hole with M — 10 M, 







Observer 



11 { 























Fig. 1. Scheme to find the relation between the emission 
time t in the local laboratory frame and the arrival time 
ia- is the point where a generic photon is emitted. L is 
the distance of P from the observer. Rt is the distance of 
the black hole from the observer, rds is the radius of the 
dyadosphere. r+ is the radius of the horizon, r (t) is the 
radius of the external surface of the fireball at time t. Rq is 

def. 

defined by i?o = -Rt — ''ds- is the angle between r [t) and 
i?T . V is a unit vector along the radial expansion velocity, 
u is a unit vector orthogonal to v oriented toward rising 
•d. k is the momentum of the photons emitted toward the 
observer. Note that we have assumed k || Rt, and this is 
certainly a very good approximation (see text, section 



So we can expand the right hand side of equation ( |l8|) in 

r(t) 
-Rt 

rit) 



powers of at first order: 



L~Rt{1 



Rt 



cos'd 



(19) 



which corresponds to neglecting the lateral displacement 
from the line of sight axis. 

Substituting (|S|) into (|l^ yields: 

rit) 



t + 



rds 



■ COS1?, 



(20) 



c c 

where we have used the fact that Rt = Rq + r^s (see 
Figure I). For r(i) we can use the following expression: 



r(t)= / V (f) dt' + rds, 
Jo 

so that equation (|2^) can be written in the form: 

-t 



(21) 



t„ = t cos tf 



rds 



(22) 

c c 
which reduces to Equation (|lj) only if v is constant and 
Tds is negligible with respect to r (t) . 

Also in Eq.(p^ we can fix t^, = t* to obtain the equa- 
tion describing the surface that emits the photons detected 
at an arrival time t*. In this case, we no longer have el- 
lipsoids of constant eccentricity ^. Since the velocity is 
strongly varying from point to point, we have more compli- 
cated surfaces like the profiles reported in Fig. ^ where at 
every point there will be a tangent ellipsoid of constant ec- 
centricity, but such an ellipsoid varies in eccentricity from 
point to point. 

For a fixed time t of emission in Eq.(|2^), the allowed 
angular interval | < cos'd < 1 leads to a corresponding 
smearing of the arrival time ta over the interval 

= 2 n^vy (23) 



6 



C.L. Bianco and R.RufRni and S.-S. Xue: Elementary spike from an EMBH Black Hole: the PEM pulse 




4e+01 1 6e+01 1 

X (cm) 



Fig. 2. Equitemporal surfaces, projected into the plane of 
Fig. |l|, with Cartesian coordinates centered on the black 
hole, for a PEM pulse arising from a black hole with 
M = 1O^M0 and Q = O.lQmax, neglecting the fireball 
thickness effects (see Fig. || and par. H), with i* going 
from 0.02620 seconds (the inner surface) to 0.02662 sec- 
onds (the outer one) with steps of 2 • 10^^ seconds. The 
observer is far away along the x axis, while y axis scale is 
expanded to see the extremely elongated surfaces, making 
the spherical surfaces centered on the black hole appear 
to be the vertical lines given in the figure. 



called the angular time scale (see e.g. Piran, 1999| ) 



5. The radiation flux with respect to arrival time 
for an infinitely thin fireball 

Eq.(|l3|), integrated with respect to gives us the value 
of F it) , the flux emitted from the PEM pulse at time 
t. Instead we want to compute F{ts), the flux detected 
at an arrival time ia- In principle, we should substitute 
Eq.(p2|) into Eq. (p^) before the integration, but we have 
no analytical relation between v (t) and t, so a numeri- 
cal integration must be performed. We used the following 
approach: 

1. We fix the value of the emission time t, and let cost? 
assume a discrete set of values between 1 and 

c 

2. For each discrete value of cosi9, we compute dF (t,??) 
and ta {t, I?). 

3. Then, we change the value of t and repeat, thus letting 
cos vary with fixed t over its allowed range, and t vary 
from the beginning of the emission to the decoupling 
time, and at each step, we compute all the values of 
dF{t,d) and U {t,i}). 

4. Now we have a large number of corresponding values of 
dF and t^,, and we have to compute F (^a) as the sum 
of all the values of dF corresponding to the same t^,- 
Unfortunately, we have computed these values numer- 
ically, and so they are not exact. Thus some dF that 
should correspond to the same tg, may correspond to 
slightly different ia, and vice versa. So we define F (ta) 



^ 6e+055 




0.0262 0.0264 0.0266 

Arrival Time (tg) (s) 



Fig. 3. Time profile for a burst from a black hole with 
M — 1O^M0 and Q = O.lQmax, neglecting fireball thick- 
ness effects. The total energy released toward the detector 
is -Etot — 3.1 • lO^^erg and the time duration (Tgo) of the 



event is Tgo ~ 7.3 • 10 
seconds (see text). 



^s. The plot is made with rc — 10 



as the sum of all the values of dF corresponding to an 
arrival time between ta and ta + '"c, where the value of 
Tc is assigned "by hand." 

5. Of course, during this integration we can see which 
values of t and give us the "same" (in the meaning 
discussed above) value of ia, and in this way we obtain 
plots like the one in Fig. ^. 

This process yields F (ti^), the so called "light curve" of 
the GRB, and this is shown in Fig. H, which reveals a very 
big problem with our computation. In fact, we see that 
the simulated time profile of the GRB seems to be "cut" 
after the peak, going to zero too fast, while experimentally 
long "tails" are observed after the peak. This is because we 
have assumed that all photons are emitted from the exter- 
nal surface of the expanding fireball, but they can also be 
emitted from a certain depth inside the surface, and this 
implies a time delay in the arrival time for photons emit- 
ted from different depths. This effect is almost negligible 
at the beginning, when the screening factor S (t) satisfies 
5 (t) <C 1, but this condition changes very rapidly during 
the expansion, when S (t) 1. The time required by a 
photon to cross the entire PEM pulse is about 10~^ sec- 
onds, when 5 (i) ~ 1, in the case reported in Fig. ^, while 
the duration of the peak we have obtained is about 10~^ 
seconds. Therefore the thickness effect cannot be neglected 
and we must change our numerical integration scheme. 

6. The radiation flux with respect to arrival time 
for a thick fireball 

In the previous section we computed the emitted fiux by 
considering the suppression factor ^ jj and assuming that 
all the radiation is emitted from the external surface of 
the expanding fireball. Now, we relax this last assump- 
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tion, and modify the estimate of the arrival time given in 
Eq. ( p2| ) by including the delay due to the thickness of the 
emission region in the fireball. We assume that, at a fixed 
emission time t, only the region of the fireball between 
r (t) and r (t) ~ X {t) is active in the photon emission, so 
the new formula is: 



V (f) dt' + r^s 



f2 



COS d - 



c 



(24) 



where < r2 < A (t) . This interval, for a fixed value of t 
and for ■(9 = 0, leads to a corresponding maximum smear- 
ing of the arrival time ta over the interval 



A 



(25) 



& 1.5e+055 
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called the intrinsic time scale. 

Due to this extra factor, the numerical integration is 
also slightly different: 

1. We pick the first value of t and divide the emitting 
part of the fireball into N emitting sub-shells, each 
of width and characterized by different values of 
the depth r2, and we assume that the flux emitted by 
each sub-shell is ^ . We pick the first sub-shell at the 
surface with r2 = and vary the value of cost? from 

to 1 with a large number of very small finite steps. 
At each step we compute dF and ia. Next we pick 
the second sub-shell, increasing r2 by one step, and 
repeat the computation of dF and varying cos-&, 
and so on for every sub-shell so that the intrinsic time 
scale of photons traveling within the emitting region is 
considered. 

2. Next we pick the second value of t and divide again 
the emitting shell of the fireball of width A (t) into N 
sub-shells of width . Note that is an increasing 
function of t, while N is fixed, so these sub-shells are 
not the same as before. We repeat the computation 
of dF and ta varying cos-!? for each sub-shell. We pick 
the third value of t and repeat, and so on until the 
decoupling time. In other words, we vary cost? at fixed 
r2 and t, and vary r2 at fixed t. 

3. The last task to perform is the same of in previous 
case, assigning to the arrival time t^ the total flux of 
all the dF corresponding to an arrival time between ta 
and ta -I- . 

The results of this new computation are shown in 
Figs. |[ H and ^, respectively regarding a black hole with 
mass equal to lO^M©, lO^M© and IOMq with charge 
Q = O.lQinax, and in Fig. ^ where the three different cases 
are shown together for comparison purposes. 

The new feature here comes from the equitemporal 
surfaces of the burst. In fact, due to the extra term r2, the 
relation between t and at fixed ta is no longer monotonic. 
Thus it can happen that photons emitted with a certain 
t, -d and r2 arrive after the ones emitted at a subsequent 
time t, at a larger angle t?, but from a lower depth r2- The 
result is that the points which in Fig. lie along curves 



Fig. 4. Time profile for a burst from a black hole with 
M = and Q = O.lQmax, considering fireball thick- 

ness effects. The total energy released towards the detector 
is -Etot — 2.4 • 10^'^erg and the time duration (Tgo) of the 
event is Tgo ~ 3.0- 10~^s. The plot is made with rc — 10^** 
seconds (see text). 
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Fig. 5. Time profile for a burst from a black hole with 
M = IO'^Mq and Q = O.lQmax, considering fireball thick- 
ness effects. The total energy released toward the detec- 
tor is Etot — 4.0 • lO'^^erg and the time duration (Tgo) 
of the event is Tgo ~ 1.1 • 10~^s. The plot is made with 
Vc — 5 ■ 10^^ seconds (see text). 

are now spread on the plane, with no possibility of finding 
clean surfaces, as can be seen in Fig. |. 

In this case, the light curves are much more similar 
to the observed ones: we can see that the first part after 
the peak shows an exponential behavior, followed by a 
tail with a power law dependence. So we suppose that the 
second part can be represented by a function like: 

F^ipi-tP-)+P3, (26) 

with pi, p2 and p^ free parameters, and the first part by: 



(27) 
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Fig. 6. Time profile for a burst from a black hole with 
M = IOMq and Q = O.lQmax, considering fireball thick- 
ness effects. The total energy released toward the detec- 
tor is ii^tot — 5.2 • lO^^erg and the time duration (Tgo) 
of the event is T90 ~ 3.9 • 10~"^s. The plot is made with 
Tc = 2 • 10~^ seconds (see text). 
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Fig. 7. The three peaks of Figs. ^, |^ and || are shown on 
the same plot, for visual comparison. 



with p4, p5 and pe free parameters. We now apply a "fit" 
algorithm separately on the two parts of the "tail" to es- 
timate the parameters. The results are shown in Figs. ^, 
1^ and 0. 

It seems that there is a relation between the values 
of p2 and p5. In fact, decreasing the value of the black 
hole mass, and so decreasing the energy of the burst, p2 
rises, making the power law decrease slower. In contrast, 
as P5 decreases the exponential part will be steeper. So the 
decrease in the energy released in the burst corresponds 
an "L" like shape, with a very steep exponential decrease, 
followed by an almost constant power-law part. 



Fig. 8. Equitemporal Surfaces for a PEM pulse from a 
black hole with M = 1O^M0 and Q — O.lQmax, consider- 
ing fireball thickness effects (see Fig. ^ and par. with 
the same scale as Fig. |2[ This plot only shows the points 
for i* = 0.03 seconds. The observer is far away along the 
X axis and the black hole is at the point (0,0). The ver- 
tical lines are a portion of spherical surfaces centered on 
the black hole. Note that the points does not form defi- 
nite lines, but show the effect of the superposition of many 
surfaces emitting independently. 

7. The temporal structure and time duration 

Within the validity of our approximation, which has the 
very strong assumption of spherical symmetry, we can 
now present the main predicted theoretical features of the 
GRBs which should be compared with observations: 

1. The "fast-rise, exponential-decay" profile which repre- 
sents a clear time trigger of the GRB signal, charac- 
terized by an exponential-decay profile with ^ e~^*'^ 
(see Figs. |, |l^ and |ri|) . 

2. The "tail" of the total radiation flux with a power law 
dependence given by ~ where a is close to unity 
(see Figs. |, |l^ and |ri|). 

From the temporal profiles of the radiation flux as a func- 
tion of the EMBH masses we flnd that the more massive 
systems have an exponential behavior with a smaller ab- 
solute value of the coefficient of the exponent and, corre- 
spondingly, a larger absolute value of the exponent of the 
power law. These effects can be compared and contrasted 
with GRB observations (see section^). 

The time duration of the GRB is usually observation- 
ally characterized by the so-called "Tgo criterion" , that is 
the time interval starting (ending) when the energy de- 
tected is the 5% (95%) of the total emitted. Namely the 
Tgo is the time duration of the emission of 90% of the en- 
ergy. This criterion is very useful in describing observed 
GRBs, because, due to the background noise, it is very 
difficult to find the exact starting time and ending time 
of the emission, while at 5% and 95% of the emission the 
signal is usually well above the noise. We can apply the 
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Fig. 9. Fit of the "tail" of the time profile of Fig. |. The Fig. 11. Fit of the "tail" of the time profile of Fig. |. The 
values obtained for the parameters are: values obtained for the parameters are: 



Parameter 


Value 


Measure Unit 


Parameter 


Value 


Measure Unit 


Pi 


(1,91 ±0,85) • lO'"' 


erg 


Pi 


(1,60 ±0,75) • 10'"'' 


erg 


s(l + I>2) 


s(l + P2) 


P2 


-1,30±0, 12 


None 


P2 


-0,900 ±0,073 


None 


P3 


(-6, 02 ±0,82) ■ 10^^ 


erg 


P3 


(-1,23 ±0,15) ■ 10'^'' 


erg 


P4 


(3, 09 ± 0, 68) • 10"** 


erg 


Pi 


(8, 5 ±2, 5) ■ 10*'" 


erg 
s 


P5 


-280,0 ±8, 3 


s-1 


P5 


(-250 ± 13) ■ 10 


s-^ 


P6 


(7, 15 ± 0, 14) • 10^" 


org 

s 


Pa 


(1,332 ±0,048) • lO^'' 


erg 
s 




Fig. 10. Fit of the "tail" of the time profile of Fig. | The 
values obtained for the parameters are: 



Parameter 


Value 


Measure Unit 


pi 


(7, 6 ± 1,8) • 10=^^ 


erg 
a(l+P2) 


P2 


-0, 999 ± 0, 046 


None 


Pa 


(-3, 16 ±0,23) ■ lO"'' 


erg 


Pi 


(1,8 ±1,1) ■10''" 


erg 
s 


P5 


-1659 ±84 


s~i 


P6 


(5, 015 ± 0,090) ■ 10'^'' 


erg 
s 



same procedure to our theoretically computed bursts, ob- 
taining their Tgo. The results for the three different black 
holes masses considered are reported in the captions of 
Figs, m P and |^. Since now the angular time scale corre- 
sponding to different emission epochs are taken into ac- 
count as well as the intrinsic time scale of the emitting 
region, the Tgo computed are significantly larger (see Fig. 

than the ones of Rufhni, et al., 1999, computed only 
on the basis of the angular time scale (see Eq.(p3|)) at the 
last transparent point. 

We emphasize the fact that by breaking the assump- 
tion of spherical symmetry that we have adopted and by 
introducing inhomogeneities in the PEM pulse, the com- 
plex sub-structures of the observed GRBs can be easily ac- 
commodated in our scenario, using the two different time 
scales (angular and intrinsic) of the emission process. An 
inhomogeneity in the PEM pulse with time scale smaller 
than the angular one (see Eq.(|23 



S < 



9 

7^c 



(i + f) 



(28) 



is invisible in the light-curve. However, any inhomo- 
geneities with time scale S within the range between the 
angular and the intrinsic time scale (see Eq.(^3|) and 
Eq.(ill)), 
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(29) 
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Fig. 12. The Tgo of the three simulated elementary bursts 
of Fig. |, I and I, plotted as a function of the black hole 
mass. 
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Fig. 13. Spectrum of the radiation observed at ta = 0.03s 
from a black hole with M = 10^ Mq and ^ = 0.1, together 
with a black body spectrum fitted on the peak. 



would give rise to complex sub-structures in the light 
curve, up to a number given by 



(30) 



8. General considerations about the spectrum 

Although the main topic of this article deals with the 
structure of the burst and its time structure as seen from 
the observer, we can add some qualitative considerations 
about the spectrum and especially on its departure from 
a black body one. If one takes the assumptions of section 
^, it is easy to compute the expected spectrum of the ob- 
served radiation. The observed number spectrum A^^, per 
photon energy e, per steradian, of photons emitted by a 
single shell is given by (in photons/eV) (see Eq.(65) of 
Ruflini, et al., 19991) : 



N,{v,T,R) EE / dV— = (5.23 x lO^^)4:nR'^dR— 
J e vj 

-ea;p[-7e(l + ^)/r]' 



log 



l-exp[-7e(l-f)/T] 



(31) 



which has a maximum at emax — 1.397r eV for 7 ^ 
1. We can then sum this spectrum over an equitemporal 
surface of our PEM-pulse to get the total spectrum of 
the radiation observed at a certain arrival time, and this 
is reported in Fig. |l^, in the case of a black hole with 
M = IO^Mq and ^ = 0.1, together with a black body 
spectrum fitted on the peak, for comparison. It is clear 
that this spectrum is already different from a black body, 
both at low and high frequency. 

These simplified results may be, however, affected by 
two additional effects, which are not very relevant for de- 
termining the equation of motion of the PEM pulse and 
the time profile of the bursts, but can indeed be of rele- 
vance for the definition of the spectrum. The first modifi- 
cation would affect some of the considerations in section 



^, where a radial dependence of the temperature of the 
radiation may be needed following a more detailed de- 
scription of the transparency condition. This will lead to 
a further broadening of the spectrum given in Fig. ^3[ 
The second modification may be due to the consideration 
of multiple inverse Compton scattering of the photons, de- 
scribed by Kompaneets equation (Kompanects, 1956; see 
also [Felten fc Rees, 1972| ). Since the increase of energy of 
a photon per scattering is: 



rUeC 



the condition for the distortion becomes 



(32) 



/kT kT 
eTiedr = 5- max (tt, t|,) > 1 (33) 
rrieC ' nieC 

with 7ie the number of the electrons, a^^e the Thomson 
cross section and tt the opacity due to Thomson scatter- 
ing 



Then, the number of scatterings is 
2 



N 



(34) 



(35) 



In the present model, at emission time, say, ~ 2 sec- 
onds after the beginning of the expansion, the comptoniza- 
tion appears to be a crucial effect. Indeed, the optical 
depth for free-free absorption is: 



rff = fc//S 10-« « 1 
where 



0.6 • 10^^ 



(36) 



(37) 
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Fig. 14. Qualitative plot of the expected spectra due to 
multiple inverse Compton scattering. 

and S is the surface density of the electrons. However, the 
medium is opaque with respect to the Compton scattering 
(described by the Thomson cross section): 



TT ~ 6 • lO'' > 1 



(38) 
i ful- 



so that the condition for comptonization (Eq.(p3|) 
filled. 

Hence, the black body spectrum, which exists in the 
comoving frame at the very initial stage of the expansion, 
quite soon will undergo modifications due to the comp- 
tonization. The observed spectrum F^,, in this case, will 
depart even more from a thermal one, with black body 
distribution at low energies, and a plateau up to the Wien 
exponential cutoff. In the Fig. |lj we give a qualitative 
form of the spectrum. 

The variation of the spectrum in the course of the ex- 
pansion of the fireball can be determined via the simulta- 
neous solution of the hydrodynamical equations of motion 



and energy balance (given in Ruffini, et al., 1999) together 
with the Kompaneets equation (Kompaneets, 1956): 



dn 
dy 



X 



1__9_ 

^ dx 



n + n -\- 



dn 
dx 



hv 
kf 



where n is the occupation number density of photons: 



(39) 



The location of the cutoff, which depends on the num- 
ber of collisions, will give independent information on the 
physical parameters of the initial Dyadosphere. This com- 
putation is beyond the scope of the present paper. 

9. Conclusions 

In our model the most gen eral GRB is c haracterized by 
six different phases (see e.g. Ruffini, 2000 ), each one with 
distinct physical processes which can lead to specific ob- 
servational features: 



The identification of the precursor. Since the mass 
range of the EMBH in which the vacuum polarization 
process can occur varies from 3.2Mq to 7.2 • IO^Mq, 
it is particularly important to identify the precursors 
of the collapsing core. While in the case of ~ IOM0 
EMBH the understanding of the process of gravita- 
tional collapse has been achieved on the ground of bi- 



nary X-ray sources (see Giacconi & Ruffini, 1978), in 
the case of more massive collapsing systems (10^ ^ 
10^ Mq) various scenarios can be considered (see e.g. 
Gurzadyan & Ruffini, in preparation). In addition to 
the identification of the precursor, it is important 
in our model to identify the magnetohydrodynamical 
conditions leading to th e charge sepa ration process in 
the collapsing core (see Ruffini, 2000 ). 
The process of gravitational collapse which leads to the 
formation of the Dyadosphere. As mentioned in section 
H, the case of spherical symmetry is currently been 
investigated ( Klippert, 2000 ). In nonspherical gravita- 
tional collapse, this phase may lead to the observa- 
tion of gravitational waves and the associated gravi- 
tationally induced electromagnetic radiation and elec- 
tromagnetically induced gravitational radiation (see 
lohnston, et al., 1973 and Johnston, et al., 1974). 
The acceleration of the sole e'^e~ electromagnetic 
plasma component (the PEM pulse) to relativistic 
Lorentz 7 factors greater than lO'^. The dynamic of 
this phase has been treated in a previous paper (see 
Ruffini, et al., 1999). This phase may last all the way 



to the reaching of the transparency condition. The ob- 
servational consequence of this phase are the topic of 
the present paper. 

The interaction of the PEM pulse, prior to reaching 
the transparency condition, with some baryonic mat- 
ter. The possible acceleration to even larger values of 
the Lorentz 7 factor (10^ ^ lO'^) of the plasma com- 
posed of 7, e"*", e~ and the electron and nucleons of the 
baryonic matter (the PEMB pulse) has been treated 
m a previous paper ( [Ruffini, et al., 20001) . The obser- 
vational aspects of this phase will be presented in a 
forthcoming paper (Bianco, Ruffini & Xue, in prepa- 
ration) where we compute the duration T90 values pre- 
dicted by this more general model as well as their tem- 
poral structure and intensity variation and compare 
and contrast these results with the observed ones. 
A substantial part of the energy of the Dyadosphere, 
in the presence of a large quantity of baryonic mat- 
ter, will be carried away by the kinetic energy of the 
baryons. The relative importance of the energy trans- 
fer from the Dyadosphere to the kinetic energy of the 
baryonic component has been estimated in a previous 



paper (Ruffini, et al., 2000). The observational aspects 



of this phase, particularly relevant to the observed af- 
terglow and possibly to neutrinos, will be presented 
in a forthcoming paper (Bianco, Fraschetti, Ruffini & 
Xue, in preparation). 

The acceleration of ultra high energy cosmic rays by 
the remnant electrodynamic structure of the EMBH 
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(see Ruffini, 199£). This problem is also under current 
examination (see Chardonnet & Ruffini, in prepara- 
tion). 

We have used the detailed model of the cre- 
ation of photons and electron-positron pairs around 
an EMBH by the process of vacuum pol ariza- 
tion (see Damour & Ruffini, 197?;, Ruffini, 1998 and 



Preparata, et al., 1998q) and t heir subsequent time evo 



lution (see Ruffini, et al., 199!: ) in order to compute the 
different relativistic effects occurring in the determination 
of the time profile of the observed radiation flux of a GRB 
with respect to arrival time, as the PEM pulse gradually 
reaches transparency. If these theoretical predictions will 
be supported by observational evidence, they will offer an 
important tool and a strong connection between the obser- 
vations of GRBs and the properties of the central engine 
which supplies the energy. This is even more compelling 
since the simplified model has only two basic parameters, 
the mass M and charge to mass ratio ^ of the EMBHs, 
giving the first clear evidence, in an astrophysical system, 
of using the extractable mass-energy of a black hole (see 



Christodoulou & Ruffini, 1971). In addition, we make pre- 



cise predictions of the structure of the bursts and of their 
time variability which can be observationally verified. 

The considerations presented in this paper refer only to 
the above mentioned phase 3 or to the initial evolution of 
phase 4. The reaching of transparency by a pure e^e~ and 
7 component can only occur for very low baryonic matter 
density around EMBH, ps ^ 10~'^g/cm?. This very spe- 
cial circumstance will define a special class of GRB with 
the specific signature presented in the present paper: short 
and elementary time variability and absence of afterglow. 
The estimate of the Tgo and the shape of the burst may 
be an important tool for their identification. The spectra 
features can also be an important tools for verifying the 
consistency of the model. 

The theoretical tools developed in this paper are ap- 
plied in the forthcoming publications to the more general 
case where baryonic matter is present. In this more general 
case only a small part of the energy of the Dyadosphere 
will be radiated in the burst. A fraction of the energy, 
increasing with the amount of baryonic matter, will be 
transferred to the kinetic energy of the baryonic matter, 
leading to observational consequences and to the afterglow 
epoch. 

Appendix A: Electron-positron pairs annihilation 
and GRBs 

The evolution of the PEM pulse is completely de- 
scribed by the general relativ istic hydrodynamical equa- 
tions (see Ruffini, et al., 1999 ) and the rate equation for 
the number-densities of electrons and positrons 



1 



(ne±[/^);^ = {n,±U'),t + ^{r'n,±U'').r 

= ^v[n^-{T)n^+{T)~n^-n^+]. 



where a is the mean annihilation-creation cross-section 
of electron-positron pairs, v is the thermal velocity of 
electron-positron pairs, 71^,+ (T) and n^- (T) are the proper 
number-densities of electrons and positrons, given by ap- 
propriate Fermi integrals. We clearly have 



{T) = n^{T), 



(A.2) 



where n^{T) is the number-density of photons given by 
the Bose integral, integrated from 2me to infinity. 

In the initial evolution stages of the PEM pulses, the 
temperature T is larger than the energy-threshold 0.5 
MeV of electron-positron pair creation (see Fig. A.l). The 
electrons and positrons created in dyadosphere are in ther- 
mal equilibrium with photons, by the process e+ + <-^- 
7 + 7- 



The rate equation becomes 
(n,±C/'^);^ = 0, 



(A.3) 



(A.4) 



that is just the conservation of the total number of 
electron-positron pairs (see Ruffini, et al., 1999 ). 

As the temperature T (see Fig. kA) drops first close to 
and then below the threshold 0.5 MeV, the electrons and 
positrons that were in thermal equilibrium with photons 
annihilate to photons. 



n^[T) > rie+c 



(A.5) 



when PEM pulses approach transparency. The ratio be- 
tween the number-densities of photons n^{T) and pairs 



n„± , defined as 



IMeV, for selected EMBHs 



, are indicated in Fig. 



for T < 



Appendix B: The screening factor in the radiation 
flux 

At the beginning, the PEM pulse is optically thick, so 
the mean free path A of the photons in the plasma (see 
Eq.(^) is more than ten orders of magnitude lower than 
the thickness D of the fireball: for a PEM pulse created 
by an EMBH with M = IO^Mq and Q = O.lQmax, the 
initial numerical values are: D ~ 10^ cm and A ~ 10~^ 
cm. So practically we have no radiation emission. During 
the expansion, the plasma density and temperature both 
decrease, while A rises, and also the emitted flux rises, 
until A ~ Z?, when the plasma is optically thin and all 
the remaining photons escape. We then assume that the 
black-body radiation flux Fbb (in erg/s) emitted by our 
flreball is given by: 



def. 



Fbb (t) ^'S (t) {t) ~ AnaT^r'cYS (t) , 



(B.l) 



where 



1.37 • 1026- 



T is the comoving frame 



(A.l) 



cni3MoV* ' 

plasma temperature given in MeV, 7 is the Lorentz 
7 factor of the bulk motion of the expanding plasma, 
r is the radius of the external surface of the PEM 
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Fig. A.l. The temperature in the comoving frame as a 
function of the local laboratory time t for typical PEM 
pulses. 

pulse and S (t) , is the "screening factor", given by (see 



Preparata, et al., 1999 ) 



S{t) = 



V2X{t) 
3 D ■ 



(B.2) 



In Fig. |B.1| are plotted the screening factors S {t) for differ- 
ent plasmas, originated around EMBHs characterized by 
charge ^ = 0.1 and by mass M = IOMq, M = IO^Mq 
and M = lO^M© respectively. 

Using the quantities obtained by the numerical 
simulation of the expansion of the PEM pulse (see 



Ruffini, ct al., 1999), we can now make a plot of the emit- 



ted flux versus emission time t, again for selected EMBHs 



(see Fig. B.2). A clear radiation flash corresponding to an 
increase of the order of magnitude of 10^° in the radiation 
flux is observed as the screening factor S (t) approaches 
unity. 

Acknowledgements. It is a pleasure to thank V. Gurzadyan, 
J. Salmonson and J. Wilson for helpful discussions, as well as 
the anonymous referee for attracting our attention on the spec- 
tral distributions, complementary to the results of the present 
work, as well as for the wording of the manuscript. 



0.01 



0.0001 



1e-006 



1e-008 



S 1e-010 



1e-012 



1e-014 



1e-016 



1e-018 



1e-020 



Mbh=10 Msoi, 
MBH=102Ms<,|a, 



Mbh=10M, 



0.0001 



0.001 




'Solar 



0.01 0.1 
t (seconds) 



10 



Fig. A. 2. The ratio between number-densities of photons 
n-y{T) and n^i, defined functions of the local 



y(r)' 



laboratory time t for the three EMBHs of Fig. B.l 



Mbh=10"Ms„|„,Q.0.1Q„ 
Mbh=102Ms„|„.Q.O.1Q„ 
Mbh=10Ms„|„.Q.0.1Q„ 



1 



10 



Time (seconds) 



Fig. B.l. The screening factor S (t) (see Eq.(p 
three different EMBHs with charge Q = O.lQmax 
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